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The  authors  wish  to  thank  Kevin  O'Neill  and  Larry  Danyluk  for  their 
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AN  AUTOMATIC  FINITE  ELEMENT  MESH  GENERATOR 

MARY  REMLEY  ALBERT  AND  JAMES  L.  WARREN 


INTRODUCTION 

The  finite  element  method  is  a  well-established  way  for  numerically 
solving  partial  differential  equations  of  Importance  in  science  and  engi¬ 
neering.  It  has  been  used  with  success  in  a  wide  variety  of  fields,  includ¬ 
ing  solid  mechanics,  heat  transfer  and  fluid  dynamics.  However,  a  draw¬ 
back  of  the  method  is  that  creating  the  initial- input  required  for  use  of 
finite  element  computer  programs  can  be  a  very  time-consuming,  tedious 
and  error-prone  procedure.  The  method  requires  the  use  of  a  mesh  that 
is  made  up  of  nodes  and  elements.  The  initial  input  to  finite  element  pro- 
grariis  requires  the  construction  of  a  list  of  all  nodes  and  their  coordinates, 
and  a  list  of  all  elements  and  nodes  that  define  the  elements,  listed  in  ,a 
particular  order  (known  as  the  incidence  list).  Haber  et  al.  (1981)  estimate 
that  roughly  80%  of  the  total  cost  of  a  finite  element  analysis  may  be  con- 
sumed^by  using. conventional  methods  for  the  input. 

In  recent  years,  the  development  of  automatic  mesh  generation  tech- 
nques  has  relieved  the  analyst  of  the  tiresome  and  costly  chore  of  creating 
the  coordinate  and  incidence  lists.  .Several  techniques  have  been  present¬ 
ed,  some  based  on  heuristic  methods,  and  others  based  on  mathematical 
concepts. 

It  is  the  purpose  of  this  paper  to  document  the  development  and  use 
of  a  user-friendly  automatic  mesh'generator  written  by  the  authors  and 
implemented  in  FORTRAN  on  the  Prime  computer  at  CRREL.  The  gener¬ 
ator  is  based  on  the  method  of  transfinite  mappings.  The  advantages  of 
the  method  will  be  discussed, .some  limitations.of  the„method -  will  be.,pre- 
sented,  and  some  notes  on  user  instruction  provided'. 


MESH  GENERATION  TECHNIQUES 

A  variety  of  mesh  generation  techniques  has  -been  documented  in  the 
literature.  Some  of  these  are  heuristic  methods,  which  may  work  in  spe¬ 
cialized  cases;  others  are  based  on  mathematical  concepts, -and  have  a 
more  predictable  outcome.  Of  the  latter  type,  three  prominent  methods 
are:  1)  Laplacian  schemes,  2)  the  use  of  isoparametric  coordinates,  and 
3)  the  use  of  transfinite  mappings.  In  this  section,. these>three  techniques 
will  be  briefly  presented. 

The  first  .technique  is -generally  termed  the  Laplacian  scheme  (Buell 
and  Bush  1973).  Essentially,  this  method  involves  the  solution  of  Laplace’s 
equation,  where  the  variables  represent  coordinate  locations  in  two-dimen- 
sionalspace: 


(1) 


a2x  a2* 
au2  +  a  v2 
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au2  av2 


o. 


(2) 


It  is  expected  that  solution  of  diffusion  equations  for  the  coordinates  will 
produce  meshes  where  the  nodes  are  distributed  evenly  on  the  interior 
of  the  grid.  Equations  1  and  2  are  then  represented  by  a  centered  finite 
difference  scheme.  The  finite  difference  equations  are: 


xi+l,j  +  xi-l,j  +  xi ,j+l  +  *i ,j-l  “  4*1  j  -  0 

ymj  +yj-u  +  3’i , j+i  +yi,j-i  -  4yi j  =  o. 


Here,  i  increments  nodes  in  the  u  direction,  and  /  increments  nodes  in  the 
v  direction.  These  equations  are  then  solved  to  find  the  coordinates  (x^j, 
yi  tj)  of  each  node  in  the  interior  of  the  grid.  This  method  generates  grids 
where  the  nodes  are  fairly  evenly  spaced  on  the  interior  of  the  grid,  but 
this  is  not  always  a  desired  effect.  There  are  variations  on  the  method. 

For  example,  Hermann  (1976)  extends  the  method  to  more  general  grid 
types,  and  combines  the  use  of  the  Laplacian  scheme  with  isoparametric 
mappings.  Denayer  (1978)  presnts  techniques  for  generating  element  con¬ 
nectivity  and  improving  the  computational  efficiency  of  the  method. 

The  other  two  methods,  isoparametric  and  transfinite,  rest  on  the  prem¬ 
ise  that  if  the  region  to  be  filled  with  a  mesh  can  be  mapped  onto  a  rec¬ 
tangular  or  triangular  space,  then  that  space  may  be  filled  with  straight 
horizontal  and  vertical  lines,  yielding  a  checkerboard  effect.  A  checker¬ 
board  requires  only  a  simple  programmed  algorithm  to  generate  a  list  of 
nodal  coordinates  and  an  .ncidence  list.  The  methods  then  map  the  checker¬ 
board  pattern  back  onto  the  original  region. 

The  use  of  isoparametric  coordinates,  as  presented  by  Zienkiewciz  and 
Phillips  (1971),  is  an  extension  of  the  use  of  isoparametric  coordinates 
for  curvilinear  l.nite  element  analysis.  Shape  functions  are  associated  with 
each  of  the  boi  uda.ry  nodes,  and  the  interior  coordinates  are  a  linear  com¬ 
bination  of  the  product  of  each  boundary  node  and  its  shape  function.  A 
major  rest:  action  on  the  method  is  that  each  side  of  the  region  must  be 
smooth;  slope  discontinuities  may  occur  only  at  the  corners. 

The  method  of  transfinite  mappings,  discussed  by  Gordon  and  Hall  (1973) 
and  Habe  et  al.  J981),  is  more  general  than  the  use  of  isoparametric  co¬ 
ordinates.  One  jf  its  greatest  advantages  is  that  slope  discontinuities  may 
occur  anywhere  along  any  boundary.  In  the  next  section,  this  method  is 
discussed  in  more  detail. 


TRANSFINITE  MAPPINGS 

This  method  of  mesh  generation  involves  the  mapping  of  a  smooth  and 
simple  region  onto  the  region  of  interest.  Because  the  mapping  guarantees 
that  the  boundary  generated  can  match  the  true  boundary  at  a  non-denum- 
erable  number  of  points,  the  method  is  termed  "transfinite." 

Consider  first  a  region  defined  by  four  sides,  where  the  sides  are  not 
necessarily  straight  or  smooth  lines,  as  illustrated  in  Figure  1.  The  object 
of  the  mapping  is  to  determine  uniquely  the  location  of  interior  nodes, 
given  the  location  of  nodes  on  the  sides  of  the  region.  To  do  this,  the  unit 
square  will  be  mapped  onto  the  four-sided  region. 
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Figure  1.  An  irregular  region  Figure  2.  Unit  square  in  (u,v)  space, 

to  be  filled  with  a  mesh. 


The  unit  square  is  described  in  a  coordinate  system  in  u  and  v  such  that 
0  <u<  1,  0  <  v  <  1,  as  illustrated  in  Figure  2.  By  assuming  that  nodes  on  op¬ 
posite  sides  of  the  square  have  the  same  u  or  v  values,  it  is  a  trivial  matter 
to  generate  coordinate  and  incidence  lists  for  the  square  by  assigning  the 
location  of  nodes  on  the  interior  of  the  square  to  points  of  intersection  of 
lines  of  constant  u  or  v.  When  the  above  process  is  complete,  this  system  of 
subdivision  may  be  mapped  to  the  irregular,  four-sided  region.  This  transfi- 
nite  mapping  is  referred  to  in  the  literature  (Gordon  and  Hall  1973.,  Haber 
et  al.  1981)  as  the  "bilinear  projector,"  and  is  defined  as  follows. 

PB(u,v)  =  (1-v)  <Pi(u)  +  vi}j  2<u)  +  (1-u)  £  i  (v)  +  u£B(v) 

-  (l-u)(l-v)  F(0,0)  -  (1-u)  vF(0,l) 

-  uvF(l,l)  -u(l-v)  F(1,0) 

where 

0  <  u  <  1,  0  <  v  <  1. 

In  eq  3,  ^(u),  ^(u),  £  i(v)  and  £g(v)  represent  the  curves  that  describe  the 
location  of  the  edges  in  x,y  space.  The  coordinates  of  the  corners  are  indi¬ 
cated  by  F(u,v).  A  mesh  generated  to  fill  the  region  depicted  in  Figure  1  is 
illustrated  in  Figure  3.  Note  that  the  mesh  on  the  interior  of  the  region  reflects 
the  boundary  shapes.  It  is  evident  that,  to  generate  a  mesh  for  use  with  finite 
elements,  tp  and  £  are  not  required  to  be  smooth  functions,  or  even  any  known 
functions.  The  mapping  requires  only  the  coordinate  values  at  the  nodal  points 
along  the  sides.  Slope  discontinuities  may  occur  anywhere  along  any  boundary. 

If  the  region  to  be  filled  with  a  mesh  has  highly  distorted  boundary  curves, 
the  interior  mesh  generated  may  lie  outside  the  boundary.  This  was  noted  by 
Gordon  and  Hall  (1973)  and  Haber  et  al.  (1981)  for  transfinite  mappings.  Note 
that  this  problem  also  occurred  in  the  use  of  isoparametric  mappings  (Hermann 
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Figure  3.  Use  of  the  bilinear  projector  to 
create  a  mesh  for  the  region  depicted  in 
Figure  1. 


a.  Region  to  be  filled  with  a 
mesh. 


b.  Overspill  with  the  chosen  corners 


I: 7&a 


1^1 


c.  Redefining  the  comers  to  eliminate 
overspill 


d.  Subdividing  the  region  to  eliminate 
overspill 


Figure  4.  Filling  a  region  with  highly  distorted  boundary  curves  with  a  mesh. 


1976).  As  an  illustration,  consider  the  region  outlined  in  Figure  4a.  When 
the  biliner  projector  is  used,  and  the  four  corners  are  as  indicated  in  Figure 
4b,  the  resultant  mesh  contains  overspill.  In  this  case,  two  remedies  are 
available.  First,  by  redefining  the  corner  locations  to  those  shown  in  Fig¬ 
ure  4c,  the  overspill  is  eliminated,  as  illustrated  in  that  figure.  The  over¬ 
spill  may  also  be  eliminated  by  dividing  the  region  into  simpler  subregions, 
which  may  then  be  filled  with  meshes  separately  and  merged  together,  as 
illustrated  in  Figure  4d. 

Haber  et  al.  (1981)  present  the  "trilinear  projector"  that  may  be  used  in 
creating  meshes  for  three-sided  regions.  Again,  slope  discontinuities  may 
occur  along  the  boundaries.  The  trilinear  projector  is  defined  as  follows. 

Pj(u,v,w)  =  \  S(v)  +  (f3^)  n(l-v)  +  (-^)n(w) 

-  wij>(o)  -  uC(o)  -vn(o)  j 
0  <  u  <  1 ,  0  <  v  <  1,  0  <  w  <  1,  u  +  v  +  w  =  1. 


The  use  of  this  projector  is  illustrated  in  Figure  5.  The  coordinates  in  the 
transformed  space,  u,  v,  and  w,  must  lie  between  0  and  1,  must  sum  to  1, 
and  run  counterclockwise  along  the  edges  around  the  three-sided  region. 

Albert  (1984)  observed  that  overspill  may  occur  also  in  regions  with 
very  regular  boundaries,  but  where  the  nodes  along  the  edges  are  highly 
bunched.  An  example  of  overspill  in  this  situation  is  depicted  in  Figure  6. 

In  Figure  7,  by  rearrangement  of  the  position  of  nodes  along  one  of  the 
edges,  we  see  that  the  situation  has  been  corrected. 

It  may  be  observed  that  the  bilinear  projector  produces  meshes  with 
four-sided  elements,  and  the  trilinear  projector  produces  meshes  with  three- 


Figure  5.  Use  of  the  trilinear  projector .  Figure  6.  Overspill  in  a  regio n  with  reg¬ 

ular  boundaries. 
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Figure  7.  Repositioning  side  nodes  to 
eliminate  overspill 


a.  Region  to  be  filled  with  a 
mesh. 


0.7  1.4  2.1  2.8  3.5  4.2  4,9 


X 

b.  Defining  the  region  as  a  four-sided 
region. 


sided  elements.  In  program  MESH, 
presented  here,  the  user  has  an  op¬ 
tion  to  have  diagonals  drawn  through 
the  four-sided  elements  generated 
from  the  bilinear  projector  to  cre¬ 
ate  triangular  elements.  Also,  in 
MESH,  the  user  has  an  option  to 
generate  eight-noded  quadrilaterals 
by  removing  the  center  node  from 
a  cluster  of  four-noded  quadrilat¬ 
erals  from  the  bilinear  projector. 

Use  of  the  trilinear  projector  will 
always  result  in  triangular  ele¬ 
ments  in  MESH. 

Now  let  us  observe  some  of  the 
properties  of  these  mappings.  First, 
the  possibility  of  allowing  slope 
discontinuities  along  any  side 
means  that  traditional  concepts 
about  the  word  "corner''  may  be 
abandoned.  For  example,  Figure 
8a  outlines  a  region  to  be  filled 
with  a  mesh.  In  Figure  8b,  the  bi¬ 
linear  projector  is  used  to  create 
the  mesh,  and  in  Figure  8c,  the 
trilinear  projector  is  used.  In 
both  cases,  the  "corners"  of  the 
region  are  marked  with  a  solid 
dot. 

Figure  9  illustrates  a  combina¬ 
tion  of  two  zones  created  using 


c.  Defining  the  region  as  a  three-sided 
region. 


Figure  8.  Use  of  both  the  bilinear  and  trilinear  projectors  to  create  a  mesh. 
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a.  Two  zones  for  the  region  b.  Use  of  the  bilinear  projector  with 

around  a  pipe.  four-noded  quadrilateral  elements  in  the 

two  zones  depicted  in  Figure  9a. 

Figure  9.  Setting  up  a  mesh  using  the  bilinear  projector  for  the  area  around  a  pipe. 


Figure  10.  Triangular  elements  for  the 
region  around  a  buried  pipe. 


Figure  11.  Eight-noded  elements  on  a 
cilrved  surface. 


the  bilinear  projector  to  set  up  a  mesh  for  a  possible  application  involving 
the  area  around  a  pipe.  In  program  MESH,  a  variety  of  options  are  made 
available  to  the  user.  Figure  10  depicts  the  use  of  triangular  elements  to 
model  the  complete  region  around  a  buried  pipe.  Figure  11  depicts  the  use 
of  eight-noded  elements.  In  Figure  12,  the  use  of  a  sine  wave  generator  for 
automatically  inputting  edge  nodes  is  displayed.  The  program  also  has  an 
option  to  allow  the  nodes  to  be  numbered,  as  depicted  in  Figure  13.  In  Fig¬ 
ure  14,  a  more  complex  figure  is  represented  by  merging  four  zones. 


Height 


Figure  12.  Use  of  the  sine  wave  option.  Figure  13.  Use  of  the  node  numbering 

option. 


Figure  14.  Example  of  merging  zones: 
a  finite  elephant. 


COMPUTER  PROGRAM,  MESH 

This  program  is  written  in  FORTRAN 
77  and  implemented  on  the  Prime  com¬ 
puter  at  CRREL.  At  the  time  of  publi¬ 
cation  of  this  report,  all  interaction  be¬ 
tween  the  analyst  and  computer  occurs 
at  the  computer  terminal.  The  program 
does  not  require  the  use  of  a  special  ter¬ 
minal,  and  a  digitizer  is  not  required. 

The  plotting  program  is  designed  for  use 
with  an  HP  7580A  plotter. 

Note  that  throughout  this  section, 
and  the  rest  of  the  report,  the  use  of 
the  word  "zone"  denotes  one  region 
that  is  filled  with  a  mesh  by  use  of 
either  the  bilinear  or  trilinear  projec¬ 
tors.  A  complicated  region  may  be 
subdivided  into  zones,  meshes  may  be 


created  for  each  zone  separately,  and  the  zones  may  then  be  joined  together, 
or  merged,  to  create  incidence  and  coordinate  lists  with  one  continuous  global 
numbering. 


MESH  is  the  main  program.  MESH  calls  subroutines  ZONE,  RATIO,  TRFIN, 
TRQUAD,  TRTRI,  RENMBR,  MERGE,  MESHPLOT  and  BANDWIRED.  It  dis¬ 
plays  the  MESH  main  menu  and  performs  virtually  all  the  tasks  that  a  user 
would  want  except  actual  mesh  construction.  From  the  main  menu  the  user 
can  1)  specify  the  number  of  zones  in  the  mesh,  2)  call  the  ZONE  subroutine, 
where  zone  construction  actually  occurs,  3)  call  the  MESHPLOT  subroutine, 
which  plots  any  previously  constructed  mesh,  4)  change  filename  for  the  next 
mesh,  5)  get  a  mesh  from  a  file  (gets  all  the  zone's  material  type,  element 
type  and  boundary  node  information),  6)  save  a  mesh  to  a  file  (saves  all  the 
zone's  material  type,  element  type  and  boundary  node  information),  7)  call 
the  subroutine  named  BANDWIRED,  which  reduces  the  bandwidth  of  a  mesh 
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and  8)  exit  the  program  and  return  to  PRIMOS  (Prime  Operating  System). 
There  are  25  subroutines  called  by  MESH;  these  subroutines  are  briefly  de¬ 
scribed  in  Appendix  A. 

The  output  of  program  MESH  consists  of  several  files,  which  contain  in¬ 
cidence  and  coordinate  lists,  plotting  information,  and  general  and  edge  in¬ 
formation  on  the  last  finite  element  mesh  that  was  created.  The  user  pro¬ 
vides  the  general  filename  during  the  run  of  the  program;  this  name  is  re¬ 
ferred  to  as  "filename"  in  the  discussion  of  output  files  that  follows. 

File  "filename. coor"  contains  the  coordinate  list  for  the  mesh.  The  first 
line  contains  a  number  that  indicates  the  number  of  nodes  in  the  entire 
mesh.  This  is  followed  by  a  list  of  the  node  numbers  and  corresponding 
coordinates. 

File  "filename. incd"  contains  the  incidence  list  for  the  mesh.  The  first 
line  contains  two  numbers  that  indicate  the  number  of  elements  in  the  en¬ 
tire  mesh,  and  the  number  of  nodes  per  element.  This  is  followed  by  a  list 
of  the  element  numbers  and  corresponding  nodes  for  each  element  (listed 
counterclockwise),  and  the  material  type  for  each  element. 

File  "filename. plot"  is  required  if  program  MESH  is  used  to  plot  the  mesh. 
The  file  contains  the  number  of  zones  in  the  mesh,  the  plotter  pen  color  for 
each  zone,  the  first  element  number  of  each  zone,  and  the  element  type 
number  for  each  zone. 

File  "filename.edge"  will  not  exist  unless  the  user  requests  that  the  infor¬ 
mation  be  saved  by  using  option  6  on  the  main  menu.  This  file  is  required  if 
a  mesh  was  only  partially  completed,  or  if  the  user  wishes  to  change  part  of 
the  mesh  at  a  later  time.  This  file  contains  the  number  of  zones  in  the  mesh, 
the  number  of  sides  in  the  zone,  the  number  of  nodes  on  sides  1  and  3  and  2 
and  4  for  each  zone,  the  material  type  number  for  each  zone,  the  element 
type  number  for  each  zone,  and  the  coordinates  of  each  corner  followed  by 
the  coordinates  of  the  nodes  on  that  side  (i.e.,  the  coordinates  of  corner  1 
with  the  coordinates  of  the  nodes  on  side  1,  etc.). 

File  "meshinfo"  contains  useful  information  for  use  in  running  program 
MESH  after  the  first  time.  It  contains  the  number  of  zones  in  the  last 
mesh  created,  the  number  1  or  2  to  indicate  whether  the  last  mesh's  edge 
information  was  or  was  not  saved,  and  the  filename  of  the  last  mesh  created. 


USER  INSTRUCTIONS 

Although  the  reader  who  wishes  to  generate  a  mesh  may  be  tempted  to 
skip  the  first  part  of  this  report,  and  begin  his  or  her  investigation  by  reading 
this  section,  it  is  our  firm  belief  that  no  computer  program  should  be  used 
completely  as  a  "black  box."  The  potential  user  should  grasp  a  basic  under¬ 
standing  of  the  method  and  brief  introduction  to  some  of  the  benefits  and 
pitfalls  before  sitting  down  at  the  computer  terminal.  This  information  is 
provided  in  the  first  four  sections  of  the  report.  Therefore,  before  using 
program  MESH,  you  are  urged  to  read  this  brief  report  in  its  entirety. 

Defining  the  boundaries  of  the  mesh  prior  to  running  MESH 

You  will  make  the  best  use  of  your  time  (and  computer  time)  if  you  make 
some  basic  decisions  before  getting  onto  the  computer.  This  planning  stage 
requires  pencil,  paper  and  some  thoughtful  reflection  on  your  part. 

First,  construct  and  label  the  coordinates  on  the  x  and  y  axes.  Sketch  the 
outline  of  the  region  to  be  filled  with  a  mesh.  If  the  region  has  more  than 
four  sides,  or  if  it  is  a  very  contorted  region,  or  if  the  region  is  composed  of 
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Figure  15.  Direction  of  numbering  of  the  nodes  on  a  side. 


more  than  one  material,  it  may  be  desirable  to  divide  the  region  into  smaller 
subregions,  called  "zones."  Each  zone  may  have  either  three  or  four  sides. 
(Note  that  for  the  mapping,  the  word  "side"  is  very  general,  and  may  include 
a  line  segment  with  a  90°  angle,  for  example.  This  concept  is  discussed  in 
earlier  sections  of  the  report.)  Number  each  zone,  starting  with  1.  Assign 
to  each  zone  a  material  type;  the  material  type  is  an  integer  that  will  be 
used  to  assign  material  properties  (such  as  density,  conductivity,  etc.)  to  the 
elements  for  the  finite  element  program. 

For  each  zone,  number  each  of  the  corners,  starting  with  1,  in  a  counter¬ 
clockwise  direction.  We  suggest  that  the  lower  left-hand  corner  of  each 
zone  be  designated  as  corner  1  for  that  zone.  For  each  zone  you  must  decide 
if  you  want  triangular  elements,  four-noded  quadrilateral  elements,  or  eight- 
noded  quadrilateral  elements. 

For  each  zone,  number  the  sides,  starting  with  1,  in  a  counterclockwise 
direction.  The  convention  for  the  numbering  is  important:  for  a  three-sided 
region,  side  1  connects  corners  1  and  3,  and  for  a  four-sided  region,  side  1 
connects  corners  1  and  4.  We  suggest  that  you  write  down  information  on 
corners  and  sides  for  reference  when  you  are  running  the  program.  Also,  the 
mesh  generation  technique  used— transfinite  mappings— relies  on  a  specific 
direction  of  numbering  for  the  nodes  on  a  side  (not  always  a  counterclock¬ 
wise  direction).  Thus,  it  will  be  important  for  you  to  be  aware  of  the  direc¬ 
tion  of  each  side,  especially  for 
sides  that  two  zones  have  in  common. 
The  convention  for  this  is  illustrated 
in  Figure  15. 

Check  to  be  sure  that  a  side  that 
is  shared  by  two  zones  will  run  in  a 
direction  consistent  with  the  conven¬ 
tion  for  each  zone.  As  an  example, 
Figure  16  depicts  the  directions  for 
each  side  for  the  mesh  illustrated  in 
Figure  9b.  Note  that  side  3  of  zone 
1  and  side  1  of  zone  2  are  shared 
sides,  yet  the  direction  is  compati¬ 
ble  with  the  convention  for  each 
zone.  (Note  that  it  is  possible  to 
match  sides  whose  nodes  run  in  dif¬ 
ferent  directions,  but  in  that  case 
the  nodes  will  have  to  be  matched 
one  by  one.) 

Determine  how  the  edge  nodes 


Figure  16.  An  example  of  a  side 
shared  by  two  zones. 
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on  each  side  of  each  zone  are  to  be  constructed,  as  follows.  First,  decide 
how  many  nodes  are  to  be  on  each  side,  remembering  that  a  three-sided  re¬ 
gion  must  have  the  same  number  of  nodes  on  each  side,  while  in  a  four-sided 
region,  alternate  sides  must  have  the  same  number  of  nodes.  Also,  if  eight- 
noded  quadrilateral  elements  are  to  be  used  in  a  four-sided  region,  there 
must  be  an  odd  number  of  nodes  on  each  side.  Next,  decide  how  the  nodes 
should  be  spaced  along  each  side.  The  following  options*  are  available  for 
specifying  the  spacing  of  the  nodes  along  each  side: 

1.  You  specify  the  nodes. 

2.  Interpolation  along  a  straight  line. 

3.  Interpolation  along  a  circle. 

4.  Interpolation  along  a  sine  wave. 

The  first  option  allows  you  to  type  in  the  coordinates  of  each  node  on  the 
side.  If  this  option  is  to  be  selected,  you  should  make  a  list  now  of  the  nodal 
coordinates,  making  sure  that  they  are  specified  in  the  proper  direction  for 
that  side,  as  discussed  above.  For  the  remaining  options,  MESH  automatic¬ 
ally  computes  the  coordinates  of  the  edge  nodes.  The  linear  interpolation 
and  sine  wave  options  allow  for  the  following  node  spacings  along  the  side: 

1.  Evenly  spaced  nodes. 

2.  Spacing  as  the  square  root  of  distance. 

3.  Linearly  increasing  spacing. 

4.  Linearly  decreasing  spacing. 

5.  A  more  severe  case  3. 

6.  A  more  severe  case  4. 

If  option  5  or  6  is  chosen,  you  are  given  a  choice  of  how  clustered  the  nodes 
are  to  be  along  the  line  by  specifying  the  fraction  of  the  line  that  is  to  be 
occupied  by  the  first  element.  The  remaining  line  segment  is  then  divided 
up  by  the  program  in  either  a  linearly  increasing  or  linearly  decreasing 
fashion,  as  indicated. 

After  you  have  determined  the  shape  of  the  mesh  and  decided  on  the  de¬ 
tails  of  the  sides  and  corners  of  each  zone,  it  is  time  to  go  to  the  computer 
and  execute  the  program  MESH.  MESH  does  not  require  the  use  of  a  special 
terminal  or  digitizer;  it  was  designed  to  run  interactively  on  any  terminal. 

The  next  section  describes  the  procedure. 

Generating  a  finite  element  mesh  using  MESH 

The  mesh-generating  program,  MESH,  is  menu-driven,  affording  the  user 
much  flexibility  in  the  procedure.  In  this  section,  each  menu  will  be  pre¬ 
sented,  with  a  brief  description  of  the  options  in  each  menu. 

Main  menu 

When  the  program  is  run,  the  first  menu  to  appear  is  the  main  menu  (shown 
on  the  following  page).  From  this  menu  you  can  initiate  the  formation  of  a 
new  mesh  (options  1,2),  retrieve  a  previously  formed  mesh  (option  5);  work 
on  a  mesh  (options  4,6,7),  plot  the  mesh  (option  3),  or  exit  to  the  operating 
system  (option  8). 

If  you  wish  to  create  a  new  mesh,  proceed  by  selecting  option  1,  followed 
by  option  2.  The  second  option  causes  the  zone  menu  to  appear,  which  allows 
the  user  to  perform  all  of  the  operations  necessary  for  the  creation  of  a 
mesh  in  each  zone.  Option  3,  "Plot  a  mesh,"  is  used  for  plotting  the  mesh 

*  Shown  exactly  as  they  are  produced  by  MESH. 
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MESH  MAIN  MENU 


1*  Specify  //  of  zones  in  the  mesh.  Currently  not  yet  defined. 

2.  Input  the  zones. 

3.  Plot  a  mesh. 

4.  Change  the  name  of  the  mesh.  Currently:  "SAMPLE 

5.  Get  a  mesh  from  a  fiLe. 

6.  Save  the  information  on  the  edges  of  the  mesh  "SAMPLE 

7.  Reduce  the  bandwidth  of  a  mesh. 

8.  Exit  to  Primos  Operation  System. 

on  the  HP  7580A  plotter.  Normally,  each  material  in  the  mesh  is  plotted  in 
a  different  color.  Also,  the  mesh  may  be  plotted  with  or  without  node  num¬ 
bers.  The  plot  option  may  be  called  either  if  a  mesh  has  just  been  created 
by  MESH,  or  if  the  incidence,  coordinate  and  plotfile  lists  for  a  mesh  al¬ 
ready  exist,  i.e.,  if  MESH  has  been  run  previously. 

The  fourth  option  in  the  main  menu,  "Change  the  name  of  the  mesh,"  is 
used  to  specify  the  names  of  the  incidence,  coordinate  and  plot,  and  edge 
information  files  of  the  mesh  to  be  generated  by  this  run  of  the  program. 

The  default  mesh  name  is  that  for  the  previous  mesh;  this  name  is  listed  in 
the  option.  The  resulting  incidence  and  coordinate  lists  will  have  the  file¬ 
name  as  specified  by  option  4  with  ".incd"  and  ".coor"  extensions  on  the 
name.  A  file  needed  by  the  plotting  subroutine  will  have  the  same  filename 
with  a  ".plot"  extension.  If  you  opt  to  save  the  information  on  edges  and 
corners  of  the  mesh,  that  information  will  be  put  into  a  file  with  filename 
(as  specified  in  option  4)  with  an  ".edge"  extension. 

Option  6,  "Save  the  information  on  the  edges  of  the  mesh,"  is  used  when  a 
mesh  has  been  input  and  you  wish  to  save  the  zone's  edge  conditions.  This 
is  useful  for  editing  the  mesh  after  it  has  been  created;  you  may  alter  parts 
of  the  existing  mesh  (change  the  shape  of  a  side  of  a  zone  or  change  node 
spacing,  for  example)  without  having  to  input  the  entire  mesh  again.  Option 
6  also  allows  you  to  save  the  information  on  the  sides  and  edges  of  a  zone 
after  each  zone  has  been  completed.  This  is  useful  if  you  are  generating 
a  multi-zone  mesh  but  are  called  away  from  the  computer  before  completing 
the  entire  mesh.  Option  5,  "Get  a  mesh  from  a  file,"  allows  the  edge  condi¬ 
tions  for  a  previously  formed  mesh  (saved  using  option  6)  to  be  read  in. 

Option  7,  "Reduce  the  bandwidth  of  a  mesh,"  examines  the  coordinate  and 
incidence  lists  for  a  mesh,  and  renumbers  the  nodes  to  reduce  the  bandwidth 
of  the  mesh.  This  option  should  be  the  very  last  operation  performed  on  a 
mesh  once  it  has  been  generated.  Then,  selecting  option  8  closes  all  files 
and  exits  the  program  to  PRIMOS. 

Zone  menu 

The  zone  menu  appears  when  you  have  selected  option  2  from  the  main 
menu,  "Input  the  zones."  The  appearance  of  this  menu  changes  slightly 
throughout  the  mesh  generation  procedure  to  provide  you  with  useful  infor¬ 
mation.  The  first  time  that  it  is  called,  the  zone  menu  appears  as  seen  on 
the  top  of  the  next  page. 

If  option  6  of  the  main  menu  ("Get  a  mesh  from  a  file")  was  used  to  read 
in  ail  of  the  zone's  edge  conditions,  then  select  option  10  of  the  zone  menu. 
Either  the  zone  menu  will  reappear,  with  the  current  zone  number  being  one 
higher  than  the  last  zone  number,  or  the  program  will  return  to  the  main 
program  if  the  last  zone  number  equals  the  number  of  zones  in  the  current 
mesh.  If  the  mesh  was  not  read  from  a  file,  you  should  begin  by  selecting 
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ZONE  MENU  Current  zone  //  1 


1.  Specify  the  #  of  sides.  The  #  of  sides  =  not  yet  defined. 

2.  Input  the  corners. 

3.  Input  the  //  of  nodes.  The  #  of  nodes  =  not  yet  defined. 

4.  Input  the  option  for  generating  the  nodes  along  any  side. 

5.  Set  the  material  of  zone.  The  material  //  =  not  yet  defined. 

6.  Set  the  element  type.  The  element  type  is  not  yet  defined. 

7.  Help  on  the  zone  construction. 

8.  Change,  insert,  delete,  or  match  nodes. 

9.  Recalculate  the  edge  nodes  with  previous  side  options. 

10.  File  Zone  //  1  and  go  to  next  the  Zone  or  the  Main  Menu. 

11.  Emergency  exit  to  the  Main  Menu. 

option  1  to  enter  the  number  of  sides  in  the  current  zone  (either  3  or  4).  The 
program  will  either  prompt  you  for  the  corners,  or  if  the  current  zone  number 
is  greater  than  one,  then  the  program  will  ask  if  the  user  wants  "automatic 
matching"  to  occur. 

Automatic  matching  means  that  an  entire  or  partial  side  of  a  previously 
entered  zone  is  mapped  onto  an  entire  side  of  the  current  zone,  thus  this 
defines  not  only  the  corners  of  the  side  but  also  the  number  of  nodes  on  the 
side  and  the  coordinates  of  all  the  nodes  on  the  side.  (If  you  do  not  wish  to 
map  up  an  entire  side  of  the  current  zone,  you  should  choose  option  8  on  the 
zone  menu  later,  after  all  of  the  corners  and  nodes  on  all  the  sides  are  de¬ 
fined.)  If  automatic  matching  is  selected,  the  program  will  ask  for  the  num¬ 
ber  of  the  zone  from  which  you  are  matching,  the  number  of  the  side  from 
which  matching  will  occur,  and  the  number  of  the  side  on  the  current  zone 
that  will  be  matched  onto.  You  also  must  input  whether  or  not  the  entire 
side  of  the  previous  zone  would  be  matched  up,  or  whether  the  side  should 
be  matched  up  between  specified  nodes  on  the  previous  zone's  side. 

Option  3  of  the  zone  menu,  "Input  the  #  of  nodes,"  should  be  selected  to 
specify  the  number  of  nodes  on  each  side  of  the  zone.  The  use  of  transfinite 
mappings  requires  that  three-sided  zones  have  the  same  number  of  nodes  on 
each  side,  while  for  four-sided  zones,  alternate  sides  must  have  the  same 
number  of  nodes.  Also,  if  eight-noded  quadrilateral  elements  are  to  be  used, 
then  there  should  be  an  odd  number  of  nodes  on  each. side. 

The  fourth  option  allows  you  to  select  the  spacing  of  nodes  along  each 
side.  It  can  be  used  to  set  the  side  option  for  the  first  time,  or  it  can  be 
used  to  change  an  existing  side  option.  The  options  offered  are: 

1.  You  specify  the  nodes. 

2.  Interpolation  along  a  straight  line. 

3.  Interpolation  along  a  cLrcle. 

4.  Interpolation  along  a  sine  wave. 

If  option  1  is  selected,  you  will  be  prompted  to  type  in  the  coordinates  of 
each  node  along  the  side.  If  interpolation  along  a  straight  line  or  sine  wave 
is  selected,  you  will  be  given  the  following  choice  of  spacings: 

1.  Evenly  spaced  nodes. 

2.  Spacing  as  the  square  root  of  distance. 

3.  Linearly  increasing  spacing. 

4.  Linearly  decreasing  spacing. 

5.  A  more  severe  case  3. 

6.  A  more  severe  case  4. 
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An  explanation  of  these  options  was  given  in  the  previous  section. 

Option  5  for  the  material  type  takes  an  integer  as  input.  This  integer 
may  be  used  later  in  the  finite  element  code  to  indicate  material  properties. 
For  MESH,  this  integer  may  correspond  to  a  pen  color  when  the  mesh  is 
plotted. 

The  sixth  option  allows  for  the  specification  of  element  type.  If  the  zone 
is  three-sided,  a  message  will  appear  saying  that  three-sided  zones  must 
have  triangular  elements.  A  four-sided  zone  may  contain  elements  that  are 
triangles,  four-noded  rectangles,  or  eight-noded  rectangles. 

Option  8  allows  you  to  change  the  coordinates  of  an  existing  side  node 
and  insert,  delete  or  match  nodes  along  a  side.  In  this  case  the  matching 
allows  you  to  match  a  partial  side  of  the  current  zone  onto  some  or  all  of 
the  nodes  specified  in  a  previous  zone.  This  includes  the  special  case  where 
two  zones  may  share  only  one  node.  If  this  option  is  selected,  the  node  will 
be  displayed  as  follows: 

Node  Menu  -  Currently  Zone  //I 

1.  Change  the  coordinates  of  nodes. 

2.  Insert  nodes. 

3.  Delete  nodes. 

4.  Match  nodes  from  other  zones  onto  nodes  of  this  zone. 

5.  Return  to  the  zone  menu. 

For  any  of  these  options,  the  program  displays  the  node  numbers  and  coordi¬ 
nates  as  they  currently  exist,  then  prompts  the  user  as  to  what  must  be  done 
to  accomplish  the  changes. 

Option  9  is  used  to  recalculate  sides  with  options  previously  entered.  This 
option  is  most  useful  when  an  entire  mesh  has  been  entered  ard  plotted,  but 
some  of  the  zones  in  the  plot  have  too  many  or  not  enough  nodes  on  a  side, 
for  example.  If  this  situation  develops,  you  should  again  select  option  2 
("Input  the  zones")  on  the  main  menu.  Then,  when  the  zone  menu  comes  up, 
select  option  3  ("Change  #  of  nodes")  and  specify  the  number  of  nodes  de¬ 
sired.  Then  select  option  9  to  recalculate  all  of  the  edge  nodes  of  all  the 
sides,  using  the  previous  side  options  but  the  new  number  of  nodes  per  side. 
Option  9  may  not  be  used  if  the  mesh  has  been  read  in  from  a  file.  When 
all  of  the  nodes  on  all  of  the  sides  have  been  specified  as  desired,  select 
option  10. 

Option  1 1  allows  you  to  leave  the  zone  menu  and  go  to  the  main  menu  to 
save  the  information  collected  so  far  on  completed  zones.  This  is  useful 
if  you  are  called  away  from  the  terminal  before  the  mesh  is  complete.  After 
selecting  option  1 1  of  the  zone  menu,  you  should  select  option  6  of  the  main 
menu  to  save  the  information. 

As  the  information  in  the  zone  menu  is  specified,  the  appearance  of  the 
menu  changes  to  indicate  the  information  that  currently  exists.  For  example, 
after  options  ’,  2,  3  and  6  have  been  selected,  the  zone  menu  may  appear 
as  shown  at  the  top  of  page  15.  The  program  has  many  checks  to  try  to  en¬ 
sure  that  the  information  supplied  by  the  user  is  complete  and  in  accordance 
with  the  theory  behind  the  mapping  used. 

When  all  of  the  zones  have  been  created,  the  program  will  ask  you  to  sup¬ 
ply  brief  information  on  merging  (joining)  the  zones.  Every  two  zones  that 
share  nodes  must  be  merged,  in  order  that  the  shared  nodes  have  only  one 
node  number,  and  global  coordinate  and  incidence  lists  can  be  generated. 

Note  that  it  is  your  responsibility  to  check  the  resulting  mesh  to  be  sure 
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ZONE  MENU  Current  zone  //  1 


1. 

2. 


3. 

4. 

5. 

6. 

7. 

8. 
9. 

10. 

11. 


Specify  the  //  of  sides.  The  #  of  sides  =  4. 

Change  the  corners. 


X  Y 


Corner  1  is  1.00000 
Corner  2  is  4.50000 
Corner  3  is  6.00000 
Corner  4  is  2.50000 


1.00000 

1.00000 

3.00000 

8.50000 


Change  the  #  of  nodes. 

There  are  6  nodes  on  side  1,  and  7  nodes  on  side  2,  currently. 
Input  the  option  for  generating  the  nodes  along  any  side. 

Set  the  material  of  zone.  The  material  it  =  not  yet  defined. 

Change  the  element  type.  The  element  type  is  4  noded  quadrilateral. 
Help  on  the  zone  construction. 

Change,  insert,  delete,  or  match  nodes. 

Recalculate  the  edge  nodes  with  previous  side  options. 

File  Zone  it  1  and  go  to  the  next  Zone  or  Main  Menu. 

Emergency  exit  to  the  Main  Menu. 


that  the  edge  nodes  of  the  zones  match.  That  is,  a  node  that  defines  part 
of  the  side  of  one  zone  should  not  fall,  for  example,  somewhere  along  the 
side  of  an  element  where  no  node  is  intended  in  the  second  zone.  Also,  the 
program  does  not  assume  that  the  complete  side  of  one  zone  matches  the 
complete  side  of  another  zone;  rather,  zones  are  merged  on  a  node-by-node 
comparison.  This  allows  a  more  general  overall  mesh  formation. 


CONCLUSIONS 

A  program  has  beemdeveloped  using  the  method  of  transfinite  mappings 
to  automatically  generate  finite  element  meshes,  given  the  location  of  nodes 
on  the  boundary  of  the  mesh.  The  program  is  easy  to  use,  and  may  generate 
meshes  for  regions  involving  complicated  shapes.  Advantages  and  disadvan¬ 
tages  of  the  method,  and  a  variety  of  examples  have  been  presented. 
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APPENDIX  A:  SUBROUTINES  CALLED  BY  MESH 


ZONE 

ZONE  is  called  by  the  main  program.  ZONE  calls  subroutines  SPEC, 

LINR,  CIRC,  CALSIN,  MAZONE,  CHINDEL  and  WHICH.  ZONE  allows  the 
user  to  input  the  information  necessary  for  generating  the  mesh  in  a  zone. 
The  zone  menu  allows  the  user  to  1)  specify  the  number  of  sides,  2)  define 
the  coordinates  of  the  corners,  3)  set  the  number  of  nodes  on  the  sides,  4) 
select  the  option  for  generating  the  nodes  along  any  given  side,  5)  set  the 
material  type  of  the  zone,  6)  set  the  element  type  for  the  zone,  7)  get  help 
on  how  to  construct  zones,  8)  change,  insert,  delete  or  match  nodes,  9)  re¬ 
calculate  the  edge  nodes  of  all  the  sides  with  the  last  used  side  options, 
and  10)  file  the  current  zone  and  go  to  the  next  zone  or,  if  the  current  zone 
is  the  last  zone,  return  to  the  main  menu.  When  defining  the  coordinates 
of  the  corners,  if  the  zone  being  constructed  is  not  the  first  zone,  then  the 
subroutine  will  ask  if  the  user  would  like  to  match  up  any  entire  sides  of 
the  current  zone  from  previously  constructed  zones.  If  any  mistakes  are 
made  the  user  can  simply  reselect  anything  desired  and  then  input  the  cor¬ 
rected  data.  The  four  options  for  defining  the  coordinates  of  the  nodes 
on  any  side  of  the  zone  are  1)  specifying  the  nodes  by  hand,  2)  linear  inter¬ 
polation  along  the  side,  3)  circular  interpolation  along  the  side,  and  4)  sine 
wave  interpolation  along  the  side.  When  the  nodes  on  all  the  sides  are  speci¬ 
fied  correctly,  ZONE  lists  the  coordinates  in  a  file  named  NODES,  which 
is  used  by  the  mesh  generation  subroutines. 

MESHPLOT 

MESHPLOT  is  called  by  the  main  program.  MESHPLOT  calls  subroutines 
PLOTS,  PLOT,  NEWPEN,  SYMBOL,  AXIS,  and  NUMBER.  All  of  the  sub¬ 
routines  called  by  MESHPLOT  are  fully  documented  (Fellers  1984).  MESH 
PLOT  allows  the  user  to  plot  a  mesh.  Given  that  the  user  has  a  mesh  that 
has  been  constructed,  the  user  can  have  it  piotted  with  node  numbers  in 
the  plot,  or  without  node  numbers.  The  user  can  have  all  the  zones  in  the 
mesh  plotted  in  black  (pen  color  1)  regardless  of  each  zone's  material  type, 
or  can  have  each  zone  plotted  in  the  color  that  corresponds  to  its  material 
type.  The  user  must  specify  the  length  of  the  x  and  y  axis.  The  subroutine 
will  scale  the  plot  to  the  lengths  desired  and  make  sure  that  all  the  data 
will  be  on  the  plot.  The  user  will  be  given  the  option  to  specify  captions 
for  the  axes,  and  will  be  given  the  option  of  specifying  a  title  for  the  plot. 

BAND  WIRED 

BANDWIRED  is  called  by  the  main  program.  BANDWIRED  calls  only 
the  subroutine  BWRSUB.  BANDWIRED  reads  in  the  incidence  and  coord¬ 
inate  lists  of  a  mesh  and  then  reduces  the  bandwidth  of  the  mesh  by  calling 
BWRSUB.  It  then  writes  the  new  incidence  and  coordinates  lists  into  the 
same  files. 

BWRSUB 

BWRSUB  is  called  by  subroutine  BANDWIRED.  It  takes  node  and  element 
arrays  and  reduces  the  bandwidth  of  the  arrays  using  Collins'  (1973)  band¬ 
width  reduction  algorithm,  as  programmed  by  Sullivan  (1985).  It  checks 
to  make  sure  that  all  the  nodes  are  assigned  to  an  element.  It  returns  the 
old  and  new  half  bandwidths. 
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RATIO 

RATIO  is  called  by  the  main  program.  It  calls  subroutine  WHICH.  RATIO 
calculates  the  fraction  of  the  total  length  from  the  corner  of  each  side  to 
each  node  on  that  side.  The  information  is  stored  in  an  array  named  RAT(,), 
and  is  written  in  a  file  named  RATIOs. 

TRFIN 

TRFIN  is  called  by  the  main  program.  TRFIN  calls  subroutine  TRIANG. 
TRFIN  performs  the  bilinear  projector  transfinite  mapping  to  produce  the 
coordinates  jf  the  nodes  on  the  interior  of  the  zone,  given  the  edge  coor¬ 
dinates.  It  also  produces  the  incidence  list  for  the  zone  by  numbering  the 
elements  in  a  row-wise  pattern.  The  elements  produced  are  four-noded 
quadrilaterals. 

TRIANG 

TRIANG  is  called  by  subroutine  TRFIN.  TRIANG  triangularizes  rectangu¬ 
lar  elements.  It  also  produces  the  incidence  list  for  the  zone  by  number¬ 
ing  the  elements  in  a  pattern. 

TRQUAD 

TRQUAD  is  called  by  the  main  program.  TRQUAD  performs  the  bilinear 
projector  transfinite  mapping  to  produce  the  coordinates  of  the  nodes  on 
the  interior  of  the  zone,  given  the  edge  coordinates.  It  also  produces  the 
incidence  list  for  the  zone  by  numbering  the  elements  in  a  row-wise  pat¬ 
tern.  The  elements  produced  are  eight-noded  quadrilaterals. 

TRTRI 

TRTRI  is  called  by  the  main  program.  TRTRI  performs  the  trilinear  pro¬ 
jector  to  produce  the  interior  node  coordinates  for  a  three-sided  zone.  By 
numbering  the  elements  in  a  row-wise  pattern,  it  also  produces  the  inci¬ 
dence  list  for  the  zone. 

RENMBR 

RENMBR  is  called  by  the  main  program.  RENMBR  renumbers  the  nodes 
and  elements  of  a  zone.  The  elements  are  numbered  starting  from  where 
the  previous  zone's  element  numbers  ended.  The  node  numbers  also  start 
from  where  the  previous  zone's  node  numbers  ended. 

MERGE 

MERGE  is  called  by  the  main  program.  MERGE  interactively  asks  the 
user  which  zones  have  nodes  in  common.  It  then  does  a  search  to  find  the 
nodes  shared  by  two  or  more  zones,  and  renumbers  the  nodes  so  that  each 
node  in  the  problem  has  only  one  node  number,  its  location  is  stored  only 
once,  and  the  coordinate  and  incidence  lists  are  changed  to  be  consistent 
with  these  changes. 

MAZONE 

MAZONE  is  called  by  subroutine  ZONE.  MAZONE  matches  up  sides  of 
zones  with  previously  entered  sides  of  earlier  zones.  It  must  map  onto  an 
entire  side  of  a  zone.  It  determines  the  corners,  node  coordinates  and  the 
number  of  nodes  on  the  side  that  is  matched  onto. 
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CKOVBR 

CKOVER  is  called  by  subroutine  MAZONE.  It  checks  to  see  if  the  user 
is  trying  to  write  over  an  existing  corner. 

CKSAME 

CKSAME  is  called  by  subroutine  MAZONE.  It  checks  all  of  the  corners 
in  a  zone  to  see  if  the  value  of  a  corner  passed  is  that  of  the  other  corner. 

CKNODES 

CKNODES  is  called  by  subroutine  MAZONE.  It  checks  to  see  whether 
the  number  of  nodes  on  the  side  being  matched  are  the  same  as  the  number 
of  nodes  on  the  other  sides  of  a  three-sided  zone  that  have  also  been  matched, 
or  the  same  as  the  number  of  nodes  on  the  opposite  side  of  a  four-sided 
zone  that  has  also  been  matched. 

DISNO 

DISNO  is  called  by  subroutine  MAZONE  and  CHINDEL.  It  displays  the 
nodes  on  a  given  side  of  a  given  zone. 

CHINDEL 

CHINDEL  is  called  by  subroutine  ZONE.  It  allows  the  user  to  change, 
insert  and  delete  individual  nodes;  it  also  allows  the  user  to  match  nodes 
from  any  side  of  any  zone  onto  a  side  of  the  zone  that  is  being  constructed 
at  the  time  that  subroutine  CHINDEL  is  called. 

SPEC 

SPEC  is  called  by  subroutine  ZONE.  It  solicits  the  edge  coordinates  to 
be  specified  interactively  by  the  user.  Because  the  order  in  which  the  points 
are  listed  is  important  (so  that  u,  v  or  w  directions  are  consistent),  SPEC 
tells  the  user  the  corner  to  which  the  first  node  specified  must  be  the  clos¬ 
est. 

LINR 

LINR  is  called  by  subroutines  ZONE  and  CALSIN.  It  determines  the  coor¬ 
dinates  of  the  nodes  on  an  edge  by  interpolating  along  a  straight  line  be¬ 
tween  the  corners.  The  user  has  the  option  to  choose  several  node  spac- 
ings:  1)  nodes  evenly  spaced,  2)  node  space  increasing  as  the  square  root 
of  the  distance,  3)  nodes  spaced  in  a  linearly  increasing  fashion,  4)  nodes 
spaced  in  a  linearly  decreasing  fashion,  5)  a  more  severe  case  of  3,  6)  a  more 
severe  case  of  4. 

CIRC 

Subroutine  ZONE  calls  this  subroutine  to  calculate  the  coordinates  of  nodes 
on  an  edge  so  that  the  nodes  are  spaced  evenly  along  the  arc  of  a  circle. 

Three  points  determine  a  circle;  this  subroutine  uses  the  corners  of  the  edge 
as  two  of  the  points,  and  asks  the  user  to  specify  any  other  point  on  the 
circle. 

QUA 

Subroutines  CIRC  and  CALSIN  call  this  subroutine.  It  calculates  the  angle 
(0  to  2PI)  between  the  positive  x-axis,  and  a  line  connecting  two  points. 
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CALSIN 

Subroutine  ZONE  calls  this  subroutine  to  calculate  the  coordinates  of 
nodes  on  an  edge  with  various  spacing  options  along  a  sine  wave.  The  sine 
wave's  wavelength,  amplitude  and  starting  point  are  needed. 

SINLINR 

This  subroutine  is  called  by  CALSIN.  It  performs  a  linear  interpolation 
between  two  points  on  a  sine  wave  with  various  spacing  options. 

WHICH 

WHICH  is  called  by  subroutines  SPEC,  LINR,  CIRC,  CALSIN,  ZONE, 
MAZONE,  and  CHINDEL.  It  is  a  short  subroutine  that,  given  the  side  number 
and  the  total  number  of  sides  in  the  problem,  returns  the  corner  numbers 
that  go  with  the  side  under  consideration. 

ROTWSIG 

ROTWSIG  is  called  by  subroutine  MESHPLOT.  It  rounds  any  real  number 
to  two  significant  digits.  It  is  used  for  the  purpose  of  rounding  the  incre¬ 
ment  on  the  axes  to  two  significant  digits. 

Function  ROUNDTOSET 

Is  called  by  subroutine  MESHPLOT.  This  function  rounds  off  any  number 
to  a  specified  magnitude.  It  is  used  to  round  off  the  x  origin  and  the  y  ori¬ 
gin  to  the  least  significant  digit  in  the  axes  increment. 
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APPENDIX  B:  EXAMPLE  RUN  USING  PROGRAM  MESH 


In  this  example,  the  items  marked  with  a  star  are  input  by  the  user, 


OK/  SEG  AMESH 

-^SAMPlI*  ?ntsr  the  Name  of  the  next  Mesh. 


User  please  note: 

R?5or*  *  documents  the  techniques  used  in 

thi^  n?nnr2mnarih0r  snd  contains  a  user  guide  for 

...jin  The  user  guide .  explains  the  terminology 

*!?ls  Program  and  provides  some  examples.  It 
is  strongly  suggested  that  you  review  at  least  the 
user  guide  section  of  the  report  before  you  preceed. 


*0 


Please  enter  0  to  cont. 


MESH  MAIN  MENU 


Currently  not  yet  defined. 


1.  Specify  A  of  zones  in  the  mesh. 

2.  Input  the  zones. 

3.  Plot  a  mesh. 

4.  Change  the  name  of  the  mesh.  Currently:  "SAMPLE  1 

5.  Get  a  mesh  from  a  file. 

6.  Save  the  information  on  the  edges  of  the  mesh  "SAMPLE 

7.  Reduce  the  bandwidth  of  a  mesh. 

Exit  to  Primes  Operating  System. 


3. 


(not 


*1 


Type  Selection  A  then  return. 


+  1 


Hog  many  zones  gill  there  be? 


MESH  MAIN  MENU 


1.  Specify  A  of  zones  in  the  mesh.  Currently  a  =  1 

2.  Input  the  zones. 

3.  Plot  a  mesh. 

4.  Change  the  name  of  the  mesh.  Currently:  SAMPLE 

5.  Get  a  mesh  from  a  file. 

6.  Save  the  information  on  the  edges  of  the  mesh  "SAMPLE  "  (not 

7.  Reduce  the  bandwidth  of  a  mesh. 

3.  Exit  to  Primos  Operating  System. 


Type  Selection  A  then  return. 


ZONE  MENU  Current  zone  A  1 


1.  Specify  the  A  of  sides.  The  A  of  sides  =  not  yet  defined. 

2.  Input  the  corners. 

3.  Input  the  A  of  nodes.  The  A  of  nodes  =  not  yet  defined. 

4.  'Input  the  option  for  generating  the  nodes  along  any  side. 

5.  Set  the  material  of  the  zone.  The  material  A  =  not  yet  defined 

6.  Set  the  element  type.  The  element  type  is  not  yet  defined. 

7.  Help  on  the  zone  construction. 

8.  Chanqe/  insert/  delete/  or  match  nodes. 

9.  Recalculate  the  edge  nodes  with  previous  side  options. 

10.  File  Zone  A  1  and  go  to  the  next  Zone  or  the  Main  Menu. 

11.  Emergency  exit  to  the  Main  Menu. 


saved) 


saved) 
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+  1 


Type  Selection  4  then  return. 


Houi  many  sides  are  there  in  this  zone 

*  4 

Specify  x  and  y  coordinates  of  corner  #  1 
*0/0 

Specify  x  and  y  coordinates  of  corner  #  2 

*  A/0 

.  Specify  x  and  y  coordinates  of  corner  #  3 
*4/1.5 

»  Specify  x  and  y  coordinates  of  corner  #  4 
"ff  0  /  2 


*5 
*  7 


Specify  the  number  of  nodes/  not  counting 
Specify  the  number  of  nodes/  not  counting 


corners/ 

corners/ 


that  you  uiant  on  side  1 
that  you  want  on  side  2 


ZONE  MENU  Current  zone  #  1 


The  #  of  sides 
Y 

0.00000 
0.00000 
1.50000 
2.00000 


1.  Specify  the  #  of  sides.  The  H  of  sides  =  4. 

2.  Change  the  corners.  X 

Corner  1  is  0.00000 

Corner  2  is  4.000O0 

Corner  3  is  4.000OO 

Corner  4  is  O.O00O0 

3.  Change  the  i  of  nodes. 

There  are  5  nodes  on  side  1/  and  7  nodes  on  side  2/  currently. 

4.  Input  the  option  for  generating  the  nodes  along  any  side. 

5.  Set  the  material  of  tne  zone.  The  material  #  =  not  yet  defined. 

6.  Set  the  element  typo.  The  element  type  is  not  yet  defined. 

7.  Help  on  the  zone  construction. 

8.  Change/  insert/  delete/  or  match  nodes. 

Recalculate  the  edge  nodes  with  previous  side  options. 
cile  Zone  #  1  and  go  to  the  next  Zone  or  the  Main  Menu. 


9. 

10.  cile  Zone  ft  1  and  gc 

11.  Emergency  exit  to  the  Main  Menu. 


*4 


Type  Selection  H  then  return. 


*1 


Please  enter  the  side  #  to  specify  a  node  option. 


The  first  node  is  closest  to  corner*  1  (  - 

The  last  node  is  closest  to  corner*  4  (  0. 


0.00000) 

2.00C00) 


Node  specification  options  for  side'  *  1 


1.  You  specify  the  nodes 

2.  Interpolation  along  a  straight  line 

3.  Interpolation  along  a  circle 

4.  Interpolation  along  a  sine  wave. 
Please  specify  the  option  number. 


Linr  spacing  options. 


1.  Evenly  spaceo  nooes. 

2.  Spacing  as  the  square  root  of  distance. 

3.  Linearly  increasing  spacing. 

4.  Linearly  decreasing  spacing. 

5.  A  more  severe  case  3. 

6.  A  more  severe  case  4. 

.  Type  selection  #  and  return. 

*  4 

Node  #  1  coordinates  x  =  0.0000  y  =  0 

Node  A  2  coordinates  x  =  0.0000  y  =  1 

Node  A  3  coordinates  x  =  0.0000  y  =  ; 

Node  4  4  coordinates  x  =  0.0000  y  =  1 

Node  A.  5  coordinates  x  -  0.0000  y  =  1 


*M 


If  all  the  sides  are  to  have  the  same  option 
y  or  Y  otherwise  enter  anything  else. 


.5714 
.0476 
.4286 
.7143 
.904  8 
enter 


Type  0  then  return  to  continue. 

*0 


ZONE  MENU  Current  zone  A  1 


1.  Specify  the  A  of  sides. 

2.  Change  the  cornars. 


Corner  1  is 
Corner  2  is 
Corner  3  is 
Corner  4  is 
3.  Change  the  4 
There  are  5 


0.00000 
4.00000 
4.00000 
0.00000 
of  nodes, 
nodes  on  side  1/  and 


The  A  of  sides  = 

0.00000 

0.00000 

1.50000 

2.00000 


4. 


7  nodes  on  side  2/  currently. 
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<t.  Input  the  option  for  generating  the  nodes  along  any  side. 

5.  Set  the  material  of  tne  zone.  The  material  #  =  not  yet  defined. 

6.  Set  the  element  type.  The  element  type  is  not  yet  defined. 

7.  Help  on  the  zone  construction. 

8.  Change/  insert/  delete/  or  match  nodes. 

9.  Recalculate  the  edge  nodes  with  previous  side  options. 

10.  File  Zone  #  1  and  go  to  the  next  Zone  or  the  Main  Menu. 

11.  Emergency  exit  to  tne  Main  Menu. 


M* 

*2 


Type  Selection  #  then  return. 

Dlease  enter  the  side  #  to  specify  a  node  option 


The  first  node  is  closest  to  corner#  1  (  0.00000/ 

The  last  node  is  closest  to  corner#  2  (  4.00000/ 

Node  specification  options  for  side  #  2 


0.00000) 

0.00000) 


1.  You  specify  the  nodes 

2.  Interpolation  along  a  straight  line 

3.  Interpolation  along  a  circle 

4.  Interpolation  along  a  sine  wave. 

Please  specify  the  option  number. 

-K  3 

Specify  x  and  y  coords  of  another  point  on  the  circle. 

#  1  coordinates  x  =  0.4732  y  =  0. 

#  2  coordinates  x  =  0.9692  y  =  0. 

#  3  coordinates  x  =  1.4807  y  =  0. 

#  4  coordinates  x  =  2.0000  y  =  0. 

#  5  coordinates  x  =  2.5193  y  =  0. 

#  6  coordinates  x  =  3.0308  y  =  0. 

#  7  coordinates  x  =  3.5268  y  -  0. 

0 


Node 
Node 
Node 
Node 
Node 
Node 
Node 
Type 
+  0 


then  return  to  continue. 


2163 

3731 

4682 

5000 

4682 

3731 

2163 


1  . 
2. 


ZONE  MENU  Current  zone  #  1 


the 


of  sides, 
corners. 


The  #  of  sides  = 
t  Y 

0.00000 
0.00000 
1.50003 
2.00000 


4. 


Change  . . 

Corner  1  is  0.00000 

Corner  2  is  4.0G000 

Corner  3  is  4.00000 

Corner  4  is  0. 00000 

3.  Change  the  #  of  nodes.  _  , 

There  are  5  nodes  on  side  1/  and  7  nodes  on  side  2/  currently 

4.  Input  the  option  for  generating  the  nodes  along  any  side. 

5.  Set  the  material  of  the  zone.  The  material  #  =  H 

6.  Set  the  element  type.  The  element  type  is  not 

7.  Help  on  the  zone  construction. 

8.  Change/  insert/  delete/  or  match  nodes. 

9.  Recalculate  the  edge  nodes  with  previous  side  options. 

10.  File  Zone  #  1  and  go  to_the  next  Zone  or  the  Main  Menu. 

11.  Emergency  exit  to  the  Main  Menu. 


not  yet  defined, 
yet  defined. 


M<> 

Mi 


Type  Selection  #  then  return. 

Please  enter  the  side  #  to  specify  a  node  option. 


The  first  node  is  closest  to  corner#  2  (  4.00000/ 

The  last  node  is  closest  to  corner#  3  C  4.00000/ 

Noae  specification  options  for  side  #  3 


0.00000) 

1.50000) 


1.  You  specify  the  nodes 

2.  Interpolation  along  a  straight  line 

3.  interpolation  along  a  circle 

4.  Interpolation  along  a  sine  wave. 
Please  specify  the  option  number. 

Linr  spacing  options. 


Node 

Node 

Node- 

Node 

Node 

■*c5y"' 


1.  Evenly  spaced  nodes. 

2.  Spacing  as  the  square  root  of  distance. 

3.  Linearly  increasing  spacing. 

4.  Linearly  decreasing  spacing. 

5.  A  more  severe  case  3. 

6.  A  more  severe  case  4. 
selection  #  and  return. 


coordinates 

coordinates 

coordinates 

coordinates 

coordinates 


then  return  to  continue. 


4.0000 

4.0000 

4.0000 

4.0000 

4.0000 


0.4286 
0.7857 
1 .0714 
1.2357 
1  .4236 
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ZONE  MENU  Current  zone  #  1 


1.  Specify  the  #  of  sides.  The  tt 

2.  Change  the  corners. 


Corner  1  is 
Corner  2  is 
Corner 
Corner 


0.00000 

4.00000 

4.00000 

0.00000 


of  sides  = 
Y 

0.00000 

0.00000 

1.50000 

2.00000 


4. 


3  is 

4  is  _ 

3.  Change  the  #  of  nodes. 

There  are  5  nodes  on  side  1/  and  7  nodes  on  side  2/  currently. 

4.  Input  the  option  for  generating  the  nodes  along  any  side. 

5.  Set  the  material  of  tne  zone.  The  material  #  -  not  yet  defined. 

6.  Set  the  element  type.  The  element  type  is  not  yet  defined. 

7.  Help  on  the  zone  construction. 

3.  Change/  insert/  delete/  or  match  nodes. 

9.  Recalculate  the  edge  nodes  with  previous  side  options. 

10.  File  Zone  #  1  and  go  to  the  next  Zone  or  the  Main  Menu. 

11.  Emergency  exit  to  the  Main  Menu. 


*4 

*4 


Type  Selection  #  then  return. 

Please  enter  the  side  It  to  specify  a  node  option. 


The  first  node  is  closest  to  corner#  4  {  0.00000/ 

The  last  node  is  closest  to  corner#  3  (  4.00000/ 

Node  specification  options  for  side  #  4 


2.00000) 

1.50000) 


*4 


1.  You  specify  the  nodes 

2.  Interpolation  along  a  straight  line 

3.  Interpolation  along  a  circle 

4.  Interpolation  along  a  sine  wave. 
Please  specify  the  option  number. 


Please  enter  the  desired  wave  length  of  the  sine  wave. 
.  <1=PI  2=2*PI  .5=PI/2  ETC.) 

*1.5 

Please  enter  the  starting  point  of  this  sin  wave. 

,  (1=PI  2=2*PI  .5=PI/2  ETC.) 

*0 

Please  enter  the  amplitude  of  the  sin  wave. 

*1* 


Sin  spacing  options. 


Type 

*1 

Node 

Node 

Node 

Node 

Node 

Node 

Node 

*Syo* 


.  Evenly  spaced  nodes. 

.  Spacing  as  the  square  root  of  distance. 
.  Linearly  increasing  spacing. 

.  Linearly  decreasing  spacing. 

.  A  more  severe  case  3. 

.  A  more  severe  case  4. 
selection  #  and  return. 


coordinates 

coordinates 

coordinates 

coordinates 

coordinates 

coordinates 

coordinates 


then  return  to  continue. 


0.5844 

1.1456 

1.6682 

2.1497 

2.6017 

3.0456 

3.5054 


2.6128 

3.0398 

3.1578 

2.9478 

2.5013 

1.9895 

1.6057 


ZONE  MENU  Current  zone  #  1 


1. 

2. 


3. 

4. 

5. 

6. 
7. 
3. 
9. 

10. 

11. 


Specify  the 
hange 


Change 
Corner  1 
Corner  2 
Corner  3 
Corner  4 


#  of  sides, 
the'  corners. 


is 

is 

is 

is 


Change  the  # 
There  are  5 


o  f 


0.00000 
4.00000 
4.00000 
0.00000 
nodes . 


The  #  of  sides 
C  Y 

0.00000 
0.00000 
1.50000 
2.00000 


4. 


------  -  nodes  on  side  1/  and  7  nodes  on  side  2/  currently 

Change  the  option  for  generating  the  nodes  along  any  side. 

Set  tne  material  of  the  zone.  The  material  #  =  not  yet  del  - 
Set  the  element  type.  The  element  type  is  not  yet  defined. 

Help  on  the  zone  construction. 

Change/  insert/  delete/  or  match  nodes. 

Recalculate  the  edge  nodes  with  previous  side  options. 

File  Zone  #  1  and  go  to  the  next  Zone  or  the  Main  Menu. 

Emergency  exit  to  the  Main  Menu. 


k  Type  Selection  #  then  return. 

*10 

You  still  must  select  #5  to  specify  the  material  type  for  this  zone 
^Type  0-  'then  return  to  continue. 
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ZONE  MENU  Currant  zona  #  1 


1. 

2. 


3. 

4. 

5. 

6. 

7. 

8. 

9. 

10. 

11 . 


Specify 

change 


the  *  of  sides, 
ange  the  corners. 


Corner  1 
Corner  2 
Corner  3 
Corner  4 
Change  the  K 
There  are  5 


is 

is 

is 

is 


The  H  of  sides  = 
K  Y 

0.00000 
0.00000 
1.50000 
2.00000 


0.00000 
4.00000 
4.00000 
0.00000 

of  nodes.  _  , 

_ _  _  nodes  on  side  1/  and  7  nodes  on  side  2/  currently. 

Change  the  option  for  generating  the  nodes  along  any  side. 

Set  the  material  of  the  zone.  The  material  D  =  not  yet  defined. 
Set  the  element  type.  The  element  type  is  not  yet  defined. 

Help  on  the  zone  construction. 

Change/  insert/  delete/  or  match  nodes. 

Recalculate  the  edge  nodes  with  previous  side  options. 

File  Zone  #  1  and  go  to  the  next  Zona  or  the  Main  Menu. 

Emergency  exit  to  the  Main  Menu. 


*5 


Type  Selection  h  then  return. 


Currently  the  material  type  for  zone 


* 


1  not  yet  defined. 


The  integer  associated  uiith  a  material  may  be  used  to  index  the 
material  types  later/  in  the  finite  element  program.  The 
main  use  of  the  material  type  in  MESH  is  to  indicate  the  color 
number  of  the  pen  to  be  used  for  plotting.  The  different 
materials  in  your  problem  should  be  numbered  from  1  to  N. 

Each  zone  must  be  assigned  a  material  type. 


*1 


Please  type  in  the  material  type  for  zone  <f  1; 


ZONE  MENU  Current  zone  H  1 


1. 

2. 


Specify  the  #  of  sides, 
nar.ge  1 


the 


The  #  of  sides  =  4. 

corners.  X  Y 

0.00000  0.00000 

4.00000  0.00000 

4.00000  1.50000 

0.00000  2.00000 

of  nodes. 

.  .  _  nodes  on  side  1/  and  7  nodes  on  side  2/  currently. 

4.  Change  the  option  for  generating  the  nodes  along  any  side. 

5.  Change  the  material  of  the  zone.  The  material  #  =  1. 

Set  the  element  type.  The  element  type  is  not  yet  defined. 

Help  on  the  zone  construction. 

Change/  insert/  delete/  or  match  nodes. 

Recalculate  the  edge  nodes  uiith  previous  side  options. 

File  Zone  $  1  and  go  to  the  next  Zone  or  the  Main  Menu. 

Emergency  exit  to  the  Main  Menu. 


Change 
Corner  1  is 
Corner  2  is 
Corner  3  is 
Corner  4  is 
3.  Change  the  # 
There  are  5 


6. 

7. 

8. 

9. 

10. 

11  . 


*6 


Type  Selection  H  then  return. 


Element  type  options  for  zone  H  1 


1.  Triangular  elements. 

2.  Four  noded  Quadrilateral  elements. 

3.  Eight  noded  Quadrilateral  elements. 
(1.35a  cpocifv  the  option  number. 

*1' . 

ZONE  MENU  Current  zone  M  1 


1.  Specify  the  H  of  sides.  The  H  of 

2.  Change  the  corners.  X 

*  '  0.00000 


Corner  1  is 
Corner  2  is 
Corner 
Corner 
3.  Change  the 
There  are 


3  is 

4  is 


sides  = 
Y 

0.00000 

0.00000 

1.50000 

2.00000 


4. 


4.00000 
4.00000 
0.00000 
It  of  nodes. 

_ .  5  nodes  on  side  1/  and  7  nodes  on  side  2/  currently. 

4.  Change  the  option  for  generating  the  nodes  along  any  side. 

5.  Change  the  material  of  the  zone.  The  material  tt  -  1. 

6.  Change  the  element  type.  The  element  type  is  triangular. 

7.  Help  on  the  zone  construction. 

8.  Change/  insert/  delete/  or  match  nodes. 

9.  Recalculate  the  edge  nodes  with  previous  side  options. 

10.  File  Zone  H  1  and  go  to  the  next  Zone  or  the  Main  Menu. 

11.  Emergency  exit  to  the  Main  Menu. 


Type  Selection  #  then  return. 
^(10 

The  half  banoisidth  of  the  mesh  is 


.  Please  enter  0  to  cont. 

*0 


8. 
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MESH  MAIN  MENU 


1.  Specify  #  of  zones  in  the  mesh.  Currently  U  =  1 

2.  Input  the  zones. 

3.  Plot  a  mesh. 

A.  Change  the  name  of  the  mesh.  Currently:  "SAMPLE  " 

5.  Get  a  mesh  from  a  file. 

o.  Save  the  information  on  the  edges  of  the  mesh  "SAMPLE  "  (not  saved) 
7.  Reduce  the . bandwidth  of  a  mesh. 

3.  Exit  to  Primos  Operating  System. 


^^Type  Selection  #  then  return. 

Are  you  assigned  to  the  HP  plotter?  If  not  you  must  exit 
to  Primos/  and  then  use  ASHP  to  assign  the  HP  plotter. 
Please  enter  n  or  N  to  exit  and  anything  else  to  continue. 


*Y 


Would  you  like  to  plot  the  current  mesh? 
(y  or  Y  =  yes  and  anything  else  =  no) 


Note  please  load  tne  plotter  with  paper  before  proceeding. 

Please  enter  "1"  if  you  want  the  nodes  numbered 
and  anything  else  otherwise. 

-Kn 

Use  ASHP  to  assign  AMLC  7 
00  YOU  WISH  TO  PLOT  THIS  TIME 

*Y 

^Do  you  want  the  mesh  titled?  Enter  Y  or  y  for  yes  anything  else  for  no. 

.Please  enter  the  label  for  the  mesh. 

-^SAMPLE 

^Oo  you  want  the  axes  labeled?  Enter  Y  or  y  for  yes  anything  else  for  no. 

.Please  enter  the  label'  for  the  x-raxis. 

*X 

.  Please  enter  the  label  for  the  y-axis. 

*Y 

Oo  you  wish  to  have  all  the  zones  plotted  in  black  (pen  color  #1)  regardless 
of  the  material  type?  (y  or  Y  =  yes  and  anything  else  = 

*Y 

Please  enter  the  number  of  inchs  to  be  on  the  X-axis. 

■fto 

Please  enter  the  number  of  inchs  to  be  on  the  Y-axis. 

▼  6 


MESH  MAIN  MENU 


1.  Specify  P  of  zones  in  the  mesh.  Currently  It  =  1 

2.  Input  the  zones, 

3.  Plot  a  mesh. 

A.  Change  the  name  of  the  mesh.  Currently:  "SAMPLE 

5.  Get  a  mesh  from  a  file. 

6.  Save  the  information  on  the  edges  of  the  mesh  "SAMPLE  "  (not  saved) 

7.  Reduce  the  bandwidth  of  a  mesh. 

8.  Exit  to  Primos  Operating  System. 


*6 


Type  Selection  A  then  return. 


MESH  MAIN  MENU 


1 

2 

3 

A 


5. 

6. 
7. 
3. 


Specify  A  of  zones  in  the  mesh.  Currently  1  =  1 

Input  tho  zones. 

Plot  a  mesh. 

Change  the  name  of  the  mesh.  Currently:  "SAMPLE 
Get  a  mesh  from  a  file. 

Save  tne  information  on  the  edges  of  the  mesh  "SAMPLE 
Reduce  the  bandwidth  of  a  mesh. 

Exit  to  Primos  Operating  System. 


"  (saved) 


*  8 


Type  Selection  A  then  return. 
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Sample 


o 

o 

•vf 


o 

V 

CD 


o 

CD 

cvi 


o 

cvj 

cvi 


>. 


o 

o 

▼I 


o 

T 


o 

cu 


o 

1 


o 

CD 


r0.40 


0.20 


0.80 


1.40 


2.00 


2.60 


3.20 


3.80 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  repro¬ 
duced  below. 


Albert,  Mary  Remley 

An  automatic  finite  element  mesh  generator  /  by  Mary  Remley  Albert 
and  James  L.  Warren.  Hanover,  N.H.:  U.S.  Army  Cold  Regions  Research 
and  Engineering  Laboratory;  Springfield,  Va.:  available  from  National 
Technical  Information  Service,  1987. 

iii,  32  p.,  illus.;  28  cm.  (CRREL  Report  87-18.) 

Bibliography:  p.  15. 

1.  Computer  code.  2.  Differential  equations.  3.  Finite  element.  4.  Fluid 
dynamics.  5.  Heat  transfer.  6.  Mathematical  methods.  7.  Solid  mechanics. 
I.  Warren,  James  L.  II.  United  States  Army.  Corps  of  Engineers.  III.  Cold 
Regions  Research  and  Engineering  Laboratory.  IV.  Series:  CRREL  Report 
87-18. 


